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Abstract 

In order to analyze data on charged pions correlation channels, 7r + (2-7r~) and 7r~ (27r + ), 
we propose new interferometry approach using the Coulomb wave function. We show 
that to describe adequately data we have to introduce new parameter describing 
the contribution of 7r~(fci)7r + (fc2) ~~ > 7r_ (k2)n + (k\) process. Using this new formula 
we analyze data on 7r + (2-7r~) and tt~(2tt + ) channels at y/s = 91 GeV by DELPHI 
Collaboration, and estimate the magnitude of this new parameter as well as the 
degree of coherence. 
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1 Introduction 

In 1995 DELPHI Collaboration reported data of the 3rd order Bose-Einstein 
Correlations (BEC) and concluded that there is a genuine the 3rd order BEC 
in 37T~ channel and the effect of the 2nd order BEC in ir + (2ir~) channel [1,2]. 
The method used in their analyses is the formulation in terms of the plane 
wave amplitude. 
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From theoretical point of view, at almost the same time, formulations for the 
2nd order BEC by means of the Coulomb wave function had been investigated 
in Refs. [3,4] 4 . Moreover, formulations for the 3rd and 4th order BEC by 
means of the Coulomb wave function have been recently proposed in Refs. [6,7]. 
In particular, the formulation of Ref. [7] contains the degree of coherence (A 1 / 2 ) 
for the magnitude of the BE exchange term. The advantage of this approach 
is that it can be directly applicable to analyses of the data on 3ir~ channel 
with the CERN-MINUIT program. 5 

In the present study we investigate whether this approach using formula de- 
rived by means of the Coulomb wave function can be extended to an unlike 
charged combination channel, 7r + (27r~). We should stress that: 

1) Through our approach, the plane wave amplitude is used as a basic calcu- 
lation. At the second step the Coulomb wave function is utilized. 

2) The interferometry effect for 7r + (27r~), i.e., the squared amplitude by the 
Coulomb wave function, contains the following neutral current: 

n-ikjir+fa) — Tr-^Tr+ffcx) . (1) 

If the magnitude of this contribution is zero, we cannot explain the data on 
7r + (27r~) channel by DELPHI Collaboration. Their data ask for the finite 
magnitude which is expressed by a new parameter (5 in this paper. 

This paper is organized in the following way: In §2, a simple model for 7r + (27r~) 
channel is investigated. In §3, a model including full amplitudes (3! = 6) is 
presented. A new formula for 7r + (27r~) channel is also derived here. In §4, 
analyses of the data on 3n~ channel are performed. The estimated parameters 
are compared with those in 7r + (27r~) channel. In the final section, concluding 
remarks are given. 



2 A simple model 

Authors of Ref. [1] stressed that data on -k + (2it~) channel are described by 
the formulas of the 2nd order BEC. We shall study whether this statement is 
correct or not using the interferometry approach formulated in terms of the 
Coulomb wave function. To analyze data on BEC effect in e + e~ — > 7r + (27r~) + 
X channel, first of all we consider the simplest diagram shown in Fig. 1. 

4 For recent presentation of various aspects of BEC in high energy physics one 
should consult a recent review [5]. 

5 In application of formulas in Ref. [6] we have to assume the parameter of the 
interaction region (R) a priori. 
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Fig. 1. A thick line denotes the positive pion. Thin lines do negative pion. The cross 
mark (x) means the exchange effect due to the Bose-Einstein statistics, whose 
magnitude is expressed by the effective degree of coherent (A 1 / 2 ). 

The plane wave amplitudes (PW A) for n + (2ir~) channel are given as 



PW A(l; H ) = _i_e i< - fc i +> ' 3:i+fc2 ' 3:2+fc3 ' a;3 - ) 

PWEA(1; H ) = J_ e l ( k [ +) -xj+kz-x-i+k^xi) 



(2a) 
(2b) 



where (+) stands for the positive pion and PWEA represents the plane wave 
exchange amplitude due to the BE statistics. 

When the Coulomb wave functions is taken into account because of charged 
pions (7r+(27r-)), the above PW A and PWEA should be replaced by A(l; H — 
-) and EA(1; + - -), 



A(l; +-- -)= ^+ 2 ( Xl , x 2 )^ 2 - k - 3 (x 2 , x 3 )^+-(x 3 , Xl ) , 
EA(l;+--)=^+-( Xl , x 3 )^~ k ;(x 3 , x 2 )^+-(x 2 , x x ) , 

where ip^.^ ls defined as 



(3a) 
(3b) 



•tiLht ( x i, x j) = r (! + iVijV^^^'^Fi-i^, 1; %{k ifl3 - kij ■ rij)\ (4) 



with r 
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±ma/kij ((+) and (— ) are the like-charge combination and unlike-charge one, 
respectively). F[a, b; x] and T(x) are the confluent hypergeometric function 
and the gamma function, respectively. The parts of the plane wave of Eq. (3) 
are given as [8,9] 



^4( 1; _| ) PWA c (3/ 2 )j\k[ +) -x 1 +k 2 -x 2 +k 3 -x x ] ^ 

EA{V H ) e ( 3 / 2 ) i [ fc 5 +) - a, l+ fc 3- a; 2 + fc2-3:3] 



(5a) 
(5b) 
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Fig. 2. Analyses of data on tt + (2tt~) and it~(2it + ) channels by means of Eq. (6). 
X 2 /n.d.f. = 3.2/29. 

where the factor (3/2) is attributed to property of the Coulomb wave function 
and 3-body problem. 

The interferometry effect for the (-k + (2tt~)) channel is calculated as 




jy(7r+27r ) 



i=l 



(1 + iQs) II / P(^)d 3 x t \Fl ~ +) + AF 2 ( 



(6) 



where p(xi) stand for the source functions of particle i. We use the Gaussian 
distribution of the radius R, p(xi) = j^Wj^ ex P — ' anc ^ 



Qs = " k 2? + S - fc 3 ) 2 + (fe 3 - fci) 2 , (7a) 
= + - -)| 2 + \EA{1- + - -)| 2 ] , (7b) 

= Re[A(l; + - -)EA*(1; +--)]. (7c) 

The parameter A 1//2 is introduced in order to estimate the strength of the BE 
effect, where the magnitude of the cross mark (x) is expressed by A 1 / 2 . Of 
course A should be less than one (1), because it can be interpreted as the 
degree of coherence in quantum optics [See also Ref. [7]]. 

Our result is given in Fig. 2 and Table 1. As seen there, the magnitude A is 
larger than 1. This suggests that we have to consider also additional diagrams 
(contributions) to Fig. 1. Therefore we should seek other possible schemes for 
the description of 7r + (27r~) channel. 
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Table 1 

Analyses of vr+(27r-) and 7r-(2vr+) BEC by DELPHI Collaboration. The system- 
atic errors for all points are assumed to be ±0.05. Small normalizations (C) are 
attributed to the long range effect (1 + 7Q3). 
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Fig. 3. Six diagrams for three charged pions, because of 3! = 6. The BE statistics is 
taken into account between \k 2 ) and 713 \ks) pions. 

3 A model including six amplitudes 



As explained in the previous section, we should consider more complex dia- 
grams than Fig. 1. In this case the following PWA and PWEA in addition to 

PWA(1; H ) and PWEA(1; H ) are necessary, because of the increased 

number of diagrams (3! = 6), 



PWA{2 


- + -) = 


i(k[ +i ■X2+k2-X3+k3-x 1 ) 
c ) 


(8a) 


PWEA(2 


- + -) = 




(8b) 


PWA(3 


--+) = 


i(k[ +i -X3+k2-xi+k3-X2) 
c ) 


(8c) 


PWEA(3 


+) = 


g»C4 + •a33+fe2-a32+fe3-a3l) 


(8d) 



These additional amplitudes described by means of the Coulomb wave function 
are given as 



A(2; - + -)=V&(*2, x 3 )ip^ 2k3 (x 3 , x^+ kl (xi, x 2 ) , 
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Fig. 4. Diagram-like interpretation of the interferometry effect in Eq. (10). 
EA(2; - + -)=1>Sfc(x 2 , »3Xt(aJ3, * 2 ) , (9b) 

A(3; - -+)=V&(*3, ^iXfel^i, ^S7(^2, x 3 ) , (9c) 
EA(3; - - +)=V&(*3, a^Vfe^, a^g^i, x 3 ) • (9d) 

The interferometry effect for ir + (2ir~) channel described by Eqs. (3) and (9) 
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Fig. 5. (a) The exchange part between tt + and ir contained in F%. (b) A possibly 
effective interpretation for the shadow region in (a) is a "p°-meson-like contribu- 
tion". 
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Fig. 6. Analyses of the data on 7r + 27r~ and tt~2it + channels by means of Eq. (10). 
X 2 /n.d.f. = 1.1/29. 



PWA 



(typical PW) e ( 3 / 2 ) i [( fc i- fc 3)- a 'i+( fc 3- fc 2)-a'2+(fc2-fci)-3:3] ^ (lid) 



In the above equations for the sake of simplicity we neglect the suffix (+). In 
actual analyses, however we should fix the charge assignment of (+) in Q 3 . 
Moreover, we have to introduce a new parameter {(5) to describe the strength 
of the shadow parts in F% + ■* and \ cf., Fig. 4. A possible interpretation 

of the role of the shadow region is a "p°-meson-like contribution" occurring 
here in order to satisfy the conservation of the neutral current, cf., Fig. 5. 

The results obtained by our new formula, i.e., Eq. (10), are given in Fig. 6 
and Table 1. 



As seen in Table 1, A becomes negative for large (5 and exceeds unity for 
small p. The possibly sets of parameters ((3, A) leading to reasonable results 
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are ({3, A) ~ (0.5, 0.2), (0.4, 0.4), (0.3, 0.6) and (0.2, 0.96). To estimate the 
strength of A we shall now analyze BEC data on 3n~ + 3n + . 



4 Analyses of 3n + 3n + BEC by Coulomb wave function 



To analyze data on 3n BEC we can use the following formula presented in 
Ref. [7] 



(1 + 7Q 3 ) ft / ti*i)<P*i [Ft~ ] + \Ft~ ] + \ 3/2 Ft~ ] ] , (12) 



N BG 
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Fig. 7. Analyses of data on 3ir + and 3ir~ channels by means of Eq. (12). 
Table 2 



Analyses of 37r~ and 37r + BEC by DELPHI Collaboration. The systematic errors 
for all points are assumed to be ±0.05. Small normalizations (C) are attributed to 
the long range effect (1 + 7Q3). 
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The result obtained by this formula is given in Fig. 7 and Table 2. Notice 
that fixed A = 1 results in large x 2 . This means that an additional parameter 
is necessary here, i.e., that we should allow the degree of coherence (A) to 
vary as well. As seen in Table 2 we have found that there is a common region 
(P, A) ~ (0.28, 0.7) for 7r + (27r-) channel and A ~ 0.7 for 3n~ + 3tt+ channel. 
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5 Concluding remarks 



We have obtained the new formula for 7r + (27r~) channel in the unlike-3rd order 
BEC, introducing the degree of coherence (A) and the effective magnitude of 
neutral current (ft). We have analyzed BEC data on 7r + (27r~) and 7r~(27r + ) 
channels in e + e~ collision at yt~s = 91 GeV, using the new formula. (Notice 
that to compare the obtained here values of R with those obtained by using the 
plane wave approach, we have to multiply it by the factor 3/2, i.e., R^ 3lv ) — > 
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0.36 fm.) 



As seen from Table 2, the following choice of parameters is possible which 
leads to good fit to data, 

(f3, A) ~ (0.28, 0.7) , 

provided that the degree of coherence (A) is almost the same in both channels. 
This finite f3 suggests that there is the genuine 3rd order contribution even in 
7r + (27r~) channel. 6 To confirm the choice of this set of parameters we need 
other data at y/s = 91 GeV as well as at different energies. Moreover, our 
formula i.e., Eq. (10) can be also applied to the same kind data from heavy- 
ion collisions. 
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